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Momentum space three-body Faddeev-like equations are used to calculate elastic, transfer and 
charge exchange reactions resulting from the scattering of deuterons on 12 C and 16 O or protons on 
13 C and 17 0; 12 C and 16 are treated as inert cores. All possible reactions are calculated in the 
framework of the same model space. Comparison with previous calculations based on approximate 
methods used in nuclear reaction theory is discussed. 
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I. INTRODUCTION 

As discussed in the review article by Austern et al. [l| 
twenty years ago, three-body models of deuteron-induced 
reactions became important since the early studies of 
stripping theory where "the internal coordinates of 
the target nucleus are ignored and the only dynamically 
active variables are the coordinates, relative to the tar- 
get nucleus, of the interacting nucleon that is captured 
by the nucleus and the spectator nucleon that goes on to 
the detector". 

The present work goes back in time, recaptures the 
three-body concept of direct nuclear reactions that 
is common to continuum discretized coupled channels 
(CDCC) calculations Q and shows the results ob- 
tained by solving Faddeev/Alt, Grassberger, and Sand- 
has (AGS) equations 0, |, @ for elastic, transfer and 
breakup reactions where three-body dynamics plays a 
dominant role. In this work we attempt to calculate 
all observables using dynamical models based on energy- 
independent or energy-dependent optical potentials for 
the nucleon-nucleus interaction [6( and realistic neutron- 
proton (np) potentials such as CD-Bonn Q ■ Some exam- 
ples are shown for reactions initiated by deuterons on 12 C 
and 16 O, as well as protons on 13 C and 17 O. Although the 
use of energy-dependent potentials in three-body calcula- 
tions is not free of theoretical problems that are discussed 
below, the results we show demonstrate the possibilities 
and the shortcomings of this model; this is, above all, the 
aim of the present paper. In addition we present the ex- 
act derivation of an alternative set of equations that may 
serve as the basis for future investigations on improving 
approximate methods in nuclear reaction theory. 

Although deuteron-nucleus three-body models, includ- 
ing striping or pick up, have already been explored in 
the past in the framework of Faddeev/AGS equations 
starting with the pioneer work of Aaron and Shanley [1] 
to the more recent calculations of Alt et al. 9], all of 
them were drastically simplified. In most cases separa- 
ble interactions were used between pairs and the correct 
treatment of the Coulomb interaction was missing. This 
situation has now changed due to the recent progress 
in the description of proton-deuteron elastic scattering 
and breakup [13, HH where the Coulomb repulsion is 



fully included using the method of screening and renor- 
malization [l2|, EH together with realistic nuclear poten- 
tials. This technical development was applied to three- 
body nuclear reactions to test the accuracy of the CDCC 
method [lij] and the convergence of the multiple scat- 
tering series in the framework of the Glauber approx- 
imation II 511 and distorted-wave impulse approximation 
(DWIA) which are standard approximations used to 
describe nuclear reaction data. 

Some of the interaction models employed in this work 
and in CDCC calculations are formally identical, but in- 
stead of solving the three-body Schrodinger equation in 
coordinate space using a representation in terms of a set 
of eigenstates pertaining to a given subsystem Hamilto- 
nian, we solve the Faddeev/AGS equations in momentum 
space and obtain numerically well converged solution of 
the three-body problem for all reactions allowed by the 
chosen interactions. In Ref. [l4| we benchmarked the 
two methods and concluded that CDCC is indeed a re- 
liable method to calculate deuteron-nucleus elastic and 
breakup cross sections, but may not provide a sufficiently 
accurate solution of the three-body problem for transfer 
and breakup in one-neutron halo nucleus scattering from 
a proton target such as n Bc+p reactions. In those cases 
the comparison of CDCC results with experimental data 
may be misleading. 

In Sec. [Til we recall the Faddeev/AGS equations, in 
Sec. Mil we present the results for three dynamical mod- 
els, and in Sec. IIVI we compare them with the results of 
standard approximations used in nuclear reaction theory. 
Conclusions are given in Sec. |Vj 



II. THE THREE-BODY EQUATIONS 

Let's consider a system of three particles (a = 1,2,3) 
with kinetic energy operator Hq, interacting by means of 
two-body potentials v a (vi = V23 in the standard odd- 
man-out notation). The full resolvent 



G(Z) = {Z - H a -^va)- 1 



and the channel resolvent 

G a {Z) = (Z-H Q -v a y x 



(1) 



(2) 
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may be related through the AGS transition operator 
Up a (Z) as 

G(Z) = 6 a0 G a (Z) + G^{Z)Up a {Z)G a {Z). (3) 

The transition operator Up a {Z) satisfies the AGS equa- 
tion m 

Go\Z) + T o{Z) G a {Z)U aa {Z), 

a 

(4) 

where the summation on a runs from one to three, 5p a = 
1 — Span Gq(Z) = (Z — Ho)^ 1 is the free resolvent, and 
T a (Z) is the two-body transition matrix (t-matrix) that 
obeys the Lippmann-Schwinger equation for pair a, 

T a (Z) = v a +v a G (Z)T a (Z). (5) 

At a given energy E in the three-body center of mass 
(cm) system the on-shell matrix elements (ipp \Up a (E + 
iO)\ip a ) calculated between the appropriate channel 
states yield all the relevant elastic, inelastic and trans- 
fer ((3, a = 1, 2, 3) as well as breakup {(3 = 0) amplitudes. 
The channel state \ip a ) for a = 1,2,3 is the eigenstate 
of the corresponding channel Hamiltonian H a = Hq + v a 
with the energy eigenvalue E made up by the bound state 
wave function for pair a times a relative plane wave be- 
tween particle a and pair a. For breakup the final state 
is a product of two plane waves corresponding to the rel- 
ative motion of three free particles. 

The AGS equations (|4|) are Faddeev-like equations 
with compact kernel and therefore suitable for numer- 
ical solution; they are consistent with the correspond- 
ing Schrodinger equation and therefore provide exact de- 
scription of the quantum three-body problem. After par- 
tial wave decomposition Eq. (0J becomes a two variable 
integral equation which we solve by standard discretiza- 
tion of momentum variables and summation of the multi- 
ple scattering series by Pade method; more details can be 
found in Refs. [13, (3- As in all numerical calculations 
convergence of results has to be tested vis-a-vis number 
of included partial waves, mesh points and Pade steps. 

To include the Coulomb interaction between two 
charged particles we use the method of screening and 
renormalization [13, [12, G3 ■ The Coulomb potential is 
screened, standard scattering theory for short-range po- 
tentials is used in the form of Eq. ^ with parametric 
dependence on the screening radius R, and the renor- 
malization procedure is applied to obtain i?-independent 
results for sufficiently large R, that correspond to the 
unscreened limit. A complete review on this subject is 
presented in Ref . [ijj] together with a number of practical 
applications. 

III. THE DYNAMICAL MODELS 

In this section we set the three-body dynamics we ap- 
ply to study all the reactions initiated by deuterons on 



12 C and 16 as well as protons on 13 C and 1 O where 
12 C and 16 are considered as inert cores. 

Although in most nuclear reaction calculations the 
deuteron wave function is generated through a Gaussian 
potential fitted to the deuteron binding energy, which is 
then used to drive the np interaction in all other partial 
waves, we use the CD-Bonn 7] potential as our realistic 
interaction for all np partial waves including the deuteron 
channel. 

For the neutron-nucleus (nA) and proton-nucleus (pA) 
interactions we use the optical potentials of Watson et 
al. @ which are based on an optical model analysis of 
nucleon scattering from lp-shcll nuclei between 10 and 50 
MeV; the nucleus A is a structureless core of mass num- 
ber A. Although core excitation may be treated in the 
present three-body models, we discard such possibility at 
this time. Therefore the relevant parameters of this op- 
tical model fit are both energy and mass dependent and 
are fitted to the existing data over the energy and mass 
range. For specific nuclei and energy, one could perhaps 
obtain a better fit but, as mentioned in the Introduction, 
our goal is to explore the possibilities of a three-body 
model that can simultaneously describe all reactions al- 
lowed by the chosen interactions and leave the fine tuning 
for an improved model study. In all calculations nuclc- 
ons are considered as spin 1/2 particles and the nuclear 
cores as spin zero particles; the spin-orbit terms of the 
optical potentials are included as well as the full opera- 
tor structure of the CD-Bonn potential for the np pair. 
The calculations include np partial waves with total an- 
gular momentum I < 3, nA partial waves with orbital 
angular momentum L < 8, and pA partial waves with 
L < 20; the total three-particle angular momentum is 
J < 35. Depending on the reaction and energy, some of 
these quantum numbers cutoffs can be safely chosen sig- 
nificantly lower, leading, nevertheless, to well converged 
results. The pA channel is more demanding than the 
nA channel due to the screened Coulomb force, where 
the screening radius R sa 10 fm for the short-range part 
of the scattering amplitude is sufficient for convergence. 
The only exception are reactions leading to a final (Ap) 
bound state where R « 15 fm and a sharper screening 
is needed. With the above choice of the calculational 
parameters we obtain well converged results for all con- 
sidered observables such that all discrepancies with the 
experimental data can be attributed solely to the short- 
comings of the interaction models that are used. 



A. Model 1 - Energy-independent optical 
potentials 

In this case we use the traditional approach based on 
energy-independent optical potentials whose parameters 
are chosen at a fixed energy. For deuteron scattering 
from nucleus A the parameters for the nA and pA po- 
tentials are taken from Ref. 0] at half the lab energy of 
the deuteron projectile. For proton scattering from the 



(An) nucleus the pA parameters are taken from Ref. [6j 
at the lab energy of the proton beam and the nA param- 
eters at zero energy, where the imaginary part of the nA 
optical potential is zero. Small adjustments to these nA 
parameters are made in order to be able to reproduce 
the experimental binding energies of the ground and ex- 
cited single particle states of the (An) nucleus while all 
Pauli forbidden bound states of the resulting potential 
are removed as described in Ref. [2(J. Original @ and 
adjusted values of these parameters are given in Table [T] 
In the present model only the modified nA parameters 
are used in given partial waves leading to the single par- 
ticle states listed in Table HT1 for 13 C and 17 0; in all other 
nA partial waves we used the original parameters @ as 
well as for the pA optical potential. While in d + A scat- 
tering the pA and nA potentials are complex, in p+ (An) 
scattering only the pA potential is complex. Although in 
both cases we are dealing with the same particles, the 
Hamiltonians are different and, therefore, in d + A we 
cannot calculate d + A — > p + (An) , but in p + (An) we 
can calculate the inverse reaction p + (An) — > d + A, or 
even p + (An) — > p + (An)*, because the nA interaction 
is real, in contrast to d + A where it is complex. 





v R (nA) 


v R (pA) 


V ao (nA) 


V so (pA) 


Ref. [6] 


60.00 


60.00 


5.5 


5.5 


N- 12 C (s) 


67.50 


66.47 






iV- 12 C (p) 


61.67 


61.50 


20.38 


20.83 


iV- 12 C (d) 


66.42 


66.42 


5.5 


5.5 


iV- 16 (s) 


61.65 


60.94 






iV- 16 (d) 


61.47 


60.89 


5.4 


5.4 



TABLE I: Original parameters of the real part of the nucleon- 
micleus optical potential [f| (first line) and those adjusted to 
the energies of bound states or resonances in given partial 
waves, all in units of MeV. The strength of the central part 
is related to v R as Vr = v r + QAZA~ iri ± 27.0(iV - Z)/A - 
0.3i5c.m.O(£c.m.) and V ao is the strength of the spin-orbit part; 
see Ref. 0] for more details. 

Results for these studies are shown by the dotted 
curves (Ml) in Figs. HH3] for d + 12 C and p + 13 C and 
Figs.HHSlfor d+ 16 andp+ 17 at different energies. As 
mentioned above the results shown by the dotted curves 
in Fig. [T] (Fig. [4j are obtained with a different Hamil- 
tonian from those in Figs. [2HS1 (Figs- EHH])- In general 
the description of the data for elastic scattering is fairly 
reasonable and within what can be expected from corre- 
sponding CDCC calculations. For the transfer reactions 
p+ 13 C -» d+ 12 C andp+ 17 -> d+ 16 shown in Fig. [3] 
and Fig. [5] respectively, one gets a reasonable agreement 
with data in the forward direction (except for a scaling 
factor), but deviations from data increase for c . m . > 30°. 

The fact that traditional three-body models of d + A 
and p + (An) scattering are inconsistent with each other 
encouraged us to study other possibilities in order to shed 
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FIG. 1: (Color online) Differential cross section divided by 
Rutherford cross section for d + 12 C elastic scattering at Ed = 
30 MeV. Predictions of Model 1 (dotted curve), Model 2 (solid 
curve), and Model 3 (dashed-dotted curve) are compared with 
the experimental data are from Ref. [2l| ]. 
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FIG. 2: (Color online) Differential cross section divided by 
Rutherford cross section for p + 13 C elastic scattering at E p = 
35 MeV. Curves as in Fig. [1] The experimental data are from 
Ref. H3. 

light on the sensitivity of results to different dynamical 
approaches. 

B. Model 2 - Energy-dependent optical potentials 

The two-body t-matrix given by Eq. ([5]) enters the 
Faddeev/AGS equation (0J for the transition operator 
Uf3 a (Z). Even if the potential is energy independent, the 
pair t-matrix has to be calculated at the two-body ener- 
gies e = E — q^/2u ai where q a is the relative momentum 
between particle a and the cm of pair a that has to be 
integrated over when solving the Faddeev/AGS equation, 
u a is the corresponding particle-pair a reduced mass, and 
E is three-body energy in the cm system. Therefore in 
three-body calculations the particles in all pairs scatter 




FIG. 3: (Color online) Differential cross section for P + 1 C — » 
d + 12 C transfer at E p = 35 MeV. Curves as in Fig. Q] The 
experimental data are from Ref. [23l ]. 



FIG. 5: (Color online) Differential cross section divided by 
Rutherford cross section for p+ 17 O elastic scattering at E p = 
35.2 MeV. Curves as in Fig. \T\ The experimental data are 
from Ref. 
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FIG. 4: (Color online) Differential cross section divided by 
Rutherford cross section for d + 16 elastic scattering at Ed = 
25.4, 36.0, and 63.2 MeV. Curves as in Fig. Q] The experi- 
mental data are from Refs. [24l. |25|. 
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FIG. 6: (Color online) Differential cross section for p+ 17 — > 
d + ls O transfer at E p = 35.2 MeV. Curves as in Fig. [T] The 
experimental data are from Ref. [24| . 



at two-body energies between E and —00. In the case 
of the CD-Bonn potential np observables are described 
with x 2 /datum ~ 1 from zero np relative energy to the 
7T production threshold. The same cannot be said about 
the nA and pA optical potentials which in the previ- 
ous model were chosen at a fixed energy. Hence they 
describe the corresponding data at that energy but not 
over the broader range that is relevant for the solution of 
the three-body Faddeev/AGS equation. 

In the present model we take the full energy depen- 
dence of the optical potential such that when nA or pA 
pairs interact at a given positive relative energy, the used 
parameters of the optical potential fit elastic nA and 
pA scattering at that energy. In addition, when the 
energy becomes negative the corresponding potentials 
become real, energy-independent and support a num- 
ber of bound states that correspond to the ground and 
excited states of the (An) and (Ap) nucleus whereas 
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the Pauli forbidden states are removed. As mentioned 
before, the parameters of the energy-dependent optical 
potentials are slightly modified to obtain the experi- 
mental binding energies at zero energy as indicated in 
Table [T] for both nA or pA potentials in given partial 
waves. In addition, the binding energy of the Pauli for- 
bidden lps/2 state in 13 C and 13 N systems is fitted to 
the 12 C neutron and proton separation energy, respec- 
tively, whereas the lpi/2 binding energy in 17 and 17 F 
systems calculated with original parameters [(J is close 
to the corresponding nucleon separation energies of 16 O. 
The resulting binding energies are given in Table |TT] for 
i3 Cj i3 Nj 17q 5 an d i7 F nuclei In the case of jV- 12 C, 

where the adjusted parameters are quite different from 
the original ones, at positive energies vr is replaced by 
v R (E c . m .) = 60.0 + (v R - 60.0) exp(- J EJ c . m ./2) and V so is 
replaced by V^ Q (£; c . m .) = 5.5 + (V so - 5.5) exp(-£ c . m ./2), 
such that the potential preserves the description of the 
N- 12 C scattering data in the desired energy regime and 
remains a continuous function of the energy. Such a re- 
placement is not needed in the case of N- 16 where the 
adjusted parameters are very close to the original ones. 

Using energy-dependent pair interactions in three- 
body calculations is by no means free of theoretical com- 
plications, such as the problem of non-orthogonality of 
three-body wave functions at different energies as a re- 
sult of the absence of a Hamiltonian theory for the scat- 
tering process. This issue can be easily understood even 
at the two-body level. If the potential is energy depen- 
dent the two-body bound states and scattering states are 
not necessarily orthogonal, much like scattering states 
corresponding to different energies. Therefore complete- 
ness relations and three-particle unitarity may be at fault 
even in the presence of real interactions. Nevertheless 
present optical model fits, in particular the one by Wat- 
son et at, are rather weak in their energy dependence, 
as can be seen by the strength of the energy dependent 
coefficients vis-a-vis the energy independent parameters; 
furthermore this energy dependence is smooth over the 
energy range of the fit except perhaps the iV- 12 C spin- 
orbit interaction in p waves near e = 0. Even in this 
case we tried a different p-wave interactions and the re- 
sults are not very different as demonstrated in the Ap- 
pendix. For this reason we believe that the problems 
of non-orthogonality of wave functions and completeness 
may be sufficiently small to allow a serious considera- 
tion of this model given its notorious advantages such as 
consistent dynamics for both d + A and p + (An) scatter- 
ing and the possibility to calculate transfer reactions to 
p + (An) and n + (Ap) final states. 

Furthermore, one should keep in mind that the en- 
ergy dependence and the imaginary part of the optical 
potential have the same origin; they arise after the elimi- 
nation of active degrees of freedom, i.e., excitations, mul- 
ticonfiguration mixing, and breakup of nucleus A, from 
the considered Hilbert space as described earlier by Fesh- 
bach [26j. However, in a three-body system this leads in 
addition to an effective energy-dependent complex three- 



body potential, or, in general, to many-body potentials 
(up to n-body) in an n-body system as formally devel- 
oped by Polyzou and Redish ^27j in the framework of ex- 
act n-body theory. Well known examples are the three- 
and four-nucleon systems described within the interac- 
tion model with energy-independent two-body potentials 
allowing for an explicit excitation of a nucleon to a A 
isobar [28], |2^| which yields effective energy-dependent 
two-nucleon and many-nucleon forces that are mutually 
consistent. In the study of three- nucleon observables it 
was found that the A-isobar effect of the two-nucleon na- 
ture is often overcompensated by the three-nucleon force 
effect. Thus, also in the three-body nuclear reactions one 
could expect a similar situation for some observables, 
that is, a partial cancellation of the effects arising due 
to the energy dependence of the two-body potential and 
due to the three-body potential if the latter would be 
included in the calculations. However, when the energy 
dependence of the two-body optical potential is intro- 
duced in the usual phenomenological way, it is not clear 
at all what should be the consistent three-body potential. 
We therefore do not attempt to include an optical three- 
body potential in the present calculations, although such 
an extension of the Faddeev/AGS framework is possible. 





1*1/2 


2*1/2 


1P3/2 


lPl/2 


W5/2 


13 c 


38.022* 


1.857 


18.722* 


4.946 


1.092 


13 N 


33.864* 




15.957* 


1.944 




17 


37.213* 


3.272 


19.267* 


16.067* 


4.143 


17 F 


32.559* 


0.105 


15.561* 


12.348* 


0.600 



TABLE II: Binding energies (MeV) of the bound states cor- 
responding to the potential parameters of Table [I] Pauli for- 
bidden bound states that are removed are marked with *. 

In Fig. fTHTTI the solid curves (M2) show the results of 
the present fully energy-dependent model for all possible 
reactions at different energies. A number of interesting 
features emerge: 

a) Elastic scattering results shown in Figs. Q~J [21 IH 
and [5] differ quite strongly from Model 1 (dotted 
curves), particularly at large angles, and become 
considerably worse when compared to data. 

b) In the low angular region (8 cm . < 30°) p + 13 C — > 
d+ 12 C (Fig.© andp+ 17 -> rf+ 16 (Fig.© re- 
sults are very similar to those obtained with Model 
1 except for a small scaling factor. 

c) Fig. [7J shows new results for the transfer reactions 
d+ 12 C — ► p+ 13 C to ground state l/2~ and excited 
states 1/2+ and 5/2+ . Again up to 6 Qm . ~ 30° the 
calculation follows the data within a small scaling 
coefficient that may be associated with a spectro- 
scopic factor. In the case of the transfer to the 
ground state, solid curves in Figs. U2 and [7J have 
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FIG. 7: (Color online) Differential cross section for d+ 12 C — ► 
p+ 13 C transfer at E d = 30 MeV. Predictions of Model 2 (solid 
curve) and Model 3 (dashed-dotted curve) are compared with 
the experimental data from Ref. [30( ] . 
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FIG. 8: (Color online) Differential cross section for p + 13 C — > 
n+ 13 N reaction at E p = 35 MeV. Dashed curve is the predic- 
tion of Model 3', other curves as in Fig. [7] The experimental 
data are from Ref. [3lj . 
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FIG. 9: (Color online) Differential cross section for d + 16 — > 
p + ir O transfer at E d = 25.4 and 36.0 MeV. Curves as in 
Fig. [7] The experimental data are from Ref. |24|. 



similar shape as expected by detailed balance tak- 
ing into account the small difference in the energies. 
The calculations also reflect the qualitative features 
of the data. 

d) Figs . 151 and [TU1 show new results for the transfer re- 
actions d+ 16 — » p+ 17 to ground state 5/2 + and 
excited state 1/2+. Again the calculations describe 
the qualitative features of the data though scaling 
factors may be needed. 

e) Figs. [5] and QT] show new results for p + 13 C — » 
n + 13 N ground state 1/2" and p + 17 O -^n + 17 F 
ground state 5/2+ and excited state 1/2+. Al- 
though in the charge exchange reactions to the 
ground state the data is not described successfully, 
it is worth noting that inp+ 17 — > n + 17 F excited 
state (1/2+) the calculations are in very reasonable 
agreement with data, except for a small scaling fac- 
tor. 

It is worth noting at this point that a good description 
of elastic data beyond small angles does not seem to be 
necessary to get the right magnitude of the transfer cross 
sections at small angles since Figs. [3] and [6] show similar 
results for two distinct models that lead to very different 
results for the elastic cross sections at large angles (see 

Figs, una a ©. 
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FIG. 10: (Color online) Differential cross section for d+ 16 — > 
p+ 17 O transfer at E d = 25.4, 36.0, and 63.2 MeV. Curves as 
in Fig. [7] The experimental data are from Ref. [24l ]. 



C. Model 3 - A "hybrid" optical potential 
approach 

Having studied these two extreme dynamical model ap- 
proaches, the energy-independent and the fully energy- 
dependent, we attempt to study a combination of the 
two. Since we want that the relevant nuclei, 13 C, 13 N, 
17 0, and 17 F have the proper low-energy spectra in order 
to describe all the relevant transfer reactions discussed 
before, we use in this case a partial-wave dependent op- 
tical potential in the following way: a) For d + A reac- 
tions in iV- 12 C (7V- 16 0) s, p, and d waves (s and d waves) 
we use the energy-dependent optical potentials of Model 
2; for p + (An) reactions the pA potential in the above 
mentioned partial waves is energy-dependent as well, but 
the nA potential is taken over from Model 1 since it is 
sufficient to bind 13 C and 17 0; b) In all other partial 
waves we use the energy-independent optical potentials 
of Model 1 with a few nuances that are explained in the 
text, depending on whether we have d + A or p + (An) 
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FIG. 11: (Color online) Differential cross section for p+ 17 — > 
n + 17 F reaction at E v = 35 MeV. Curves as in Fig. |8] The 
experimental data are from Ref. [3^ ]. 



scattering. 

Since Model 1 is more absorptive than Model 2 due 
to the large impact of the imaginary part of the optical 
interactions on the elastic cross sections we expect this 
hybrid model to improve the description of the elastic 
data. 

For d + A scattering, results are shown by the dash- 
dotted curves (M3) in Figs. [TJ H [fl El and and [TO] In 
d + 12 C (d + 16 0) both nA and pA optical potentials 
are, like in Model 2, energy dependent in s, p, and d 
waves (s and d waves) while in all other partial waves 
they are energy independent with the parameters chosen 
at half the deuteron lab energy like in Model 1. The 
dash-dotted curves show a remarkable improvement vis- 
a-vis the fully energy-dependent calculations (solid lines 
in Model 2), particularly at large angles. This effect is 
visible not only in elastic scattering (Figs. Q] and [4]), but 
also in the transfer reactions d + A — > p + (An) shown 
in Figs. [TlIU and [TTJ1 where in some specific cases such as 
in Figs . 151 and [TTJ1 one gets quite reasonable description of 
the data. 

For p + (An) scattering results are again shown in 
Figs. iUll and [11] by the dash-dotted curves 
(M3). In p + 13 C (p + 17 O) the pA optical potentials are, 
like in Model 2, energy dependent in s, p, and d waves (s 
and d waves) and, in all other partial waves, are energy 
independent with the parameters chosen at the proton 
lab energy, like in Model 1. As for the nA optical poten- 
tial it is chosen as in Model 1 where, in all partial waves 
the potential is real and supports a number of single par- 
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ticle states as mentioned before. As in d + A reactions, 
we notice an improvement in the description of p + (An) 
elastic (Figs. [2] and [5|) as well as p + (An) — » d + A trans- 
fer (Figs. [3] and [75]) observables. Nevertheless, at small 
angles (9 c . m . < 30°), the differences between Models 1, 
2, and 3 are quite small indicating that the extracted 
spectroscopic factors would be of similar size as well. 

Finally for the charge transfer reactions shown in 
Figs. [8J and [TI] we add a new calculation shown by the 
dashed curves (M3') where the nA optical potential is 
energy-dependent in the s, p, and d partial waves for 13 C 
and s and d waves for 17 O like in Model 2 but is, in the 
other partial waves, energy-independent with the param- 
eters chosen according to the lab energy of the neutron 
in the inverse reaction n + ( Ap) — > p + (An) . The three 
curves shown in Figs . [51 and 1111 are not very different aside 
a scaling factor. 



IV. DISTORTED WAVE EQUATIONS 

In order to relate our calculations to the standard ap- 
proaches of nuclear reaction theory we derive an alterna- 
tive set of scattering equations. Introducing an effective 
interaction V a acting between particle a and the cm of 
pair a as shown in Fig. [T3] one may define a new resolvent 

G a (Z) = (Z-H -v a -V a )- 1 , (6) 

such that 

G a (Z) = G a (Z)+G a (Z)f a (Z)G a (Z), (7) 

f a (Z) =V a + V a G a (Z) f a (Z). (8) 

Likewise one may define a distorted wave in channel a as 

\^ a ) = (l+G a (Z)f a (Z))\^ a ). (9) 

Using the identity 

G(Z) = Gp(Z) + G (Z)[Gf(Z) - G-\Z)]G(Z), (10) 

together with Eqs. (fT|) and (JFJJ one gets 

G(Z) = Gp(Z) + Gp(Z)npG(Z), (11) 

where 

Hp = ^2hpv tr -Vp. (12) 



o 



a 



FIG. 12: V a optical interaction in channel a 



A new operator Up a (Z) relating G(Z) to G a (Z) instead 
of G a (Z), i.e., 

G(Z) = Sp a G a (Z) + Gp(Z) Up a (Z) G a (Z), (13) 

satisfies an equation 

Up a (Z) = ~8p a G- 1 (Z) + n, 3 + n p G a (Z)U aa (Z). (14) 

Its relation to the standard Faddeev/AGS operator 
Up a (Z) is obtained using Eq. ([7]) in (TT3"]) and compar- 
ing back with ([3]) as 



Up a (Z) = S 0a f a (Z) + 1 + fp(Z) G (Z) 



xUp a (Z) l + G a (Z)f a (Z) 



(15) 



which for on-shell elements reads 



(^f3\U/3a(Z)\tp a ) = 8i3 a (lpi3\f a (Z)\ll} a ) 



$p\Ue a (z)$ a ). 

(16) 



In Eq. (fT4"|) the term 5^ ct G~ 1 (Z) is zero on-shell and will 
be omitted in the following considerations. 

Using the Born approximation Up a (Z) ~ Hp for (3 ^ 
a, one gets 



(MU Pa (Z)\^ a ) ~ (faSlfiWa) 



(17) 



which corresponds to the usual distorted-wave Born ap- 
proximation (DWBA) for the transfer reactions in the 
post form. 

On the other hand, using V\ = in the case of particle 
1 colliding with pair (23) one gets f a (Z) = 0, Uu(Z) = 
Un(Z), and Q,\ — V2 + i>3 leading to 

Un(Z) = (v 2 + v 3 ) + (v 2 + v 3 ) Gx(Z) U U (Z), (18) 

which is the integral form of the CDCC differential equa- 
tion. This equation by itself is not connected in all orders 
of iteration and therefore cannot be solved by standard 
numerical methods since it does not satisfy the Fredholm 
alternative. Nevertheless one may follow the momentum 
space version of the CDCC approach and use the spec- 
tral decomposition of G\(Z) to obtain a set of coupled 
equations involving the continuum wave functions of pair 
(23) in addition to the bound state wave function If 
the continuum is discretized and the corresponding wave 
functions normalized a la CDCC, the solution of Eq. (TTS")) 
includes the bound to continuum and continuum to con- 
tinuum couplings that are common to CDCC calcula- 
tions. In Ref. [14| we have shown that CDCC calcula- 
tions for deuteron elastic scattering and breakup from a 
heavier target are reliable, but transfer and breakup re- 
actions involving the scattering of a halo nucleus from a 
light target, such as 11 Be + p, may be at fault. Therefore 
we expect all deuteron elastic scattering results shown 
in Figs. Q]and|l]to agree well with those obtained from 
equivalent CDCC calculations. 
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The present derivations may be useful in future stud- 
ies of approximate methods often used in nuclear reaction 
calculations involving deuterons or halo nuclei. Since at 
this time in the present framework we do not have the 
means to test the validity of Eq. (fT7)) or any other ap- 
proximation, we compare our results with published cal- 
culations involving either DWBA, coupled-channel Born 
approximation (CCBA), or various adiabatic approaches 
[H, HH HE] • These calculations use wave functions \ip a ) 
and optical potentials V a that are tuned at the consid- 
ered reaction energies, while our calculations use global 
fits to nuclear reaction data and are aimed at providing 
a description of the data in different channels simultane- 
ously. That tuning may be, at least in part, the reason 
for a better description of the data as discussed below. 

The data points in Fig. [3] for p + 13 C -> d + 12 C were 
analyzed in Ref. [23| using DWBA, adiabatic deuteron 
breakup approximation (ADBA), and CDCC-CCBA. Up 
to 30° our results coincide with all of the previous calcu- 
lations, but at larger angles adiabatic and CDCC-CCBA 
calculations follow the data much better than ours. 

Likewise, the data points in Fig.[7]for d+ 12 C — > p+ 13 C 
were analyzed in Ref. [3(| using DWBA for transition to 
the 13 C ground state and CCBA for transitions to the 
1/2+ and 5/2+ excited states. While DWBA provides a 
better fit to the data at small angles, it still overshoots 
the data at larger angles much like our calculations. As 
for the reactions leading to the l/2 + and 5/2 + excited 
states of 13 C, both CCBA and our calculations describe 
the data equally badly. 

In Ref. [31| the p+ 13 C — > n+ 13 N reaction was analyzed 
with DWBA which undershoots the data at small angles, 
much like our results shown in Fig.[8l but, overall, provide 
a better description of the data. 

The results in Fig. [Hand Fig. QUI for d+ 16 -> p + 17 O 
may be compared with DWBA studies from Ref. [2J] and 
adiabatic calculations from Ref. [35| for 36 MeV and 63.2 
MeV deuterons. The results we get in the framework 
of Model 3 (the "hybricT) are qualitatively similar to 
those obtained in Ref. 0, [35| although quantitatively 
they may differ in specific angular regions, leading to a 
description of the data that is not as good as the one 
provided by DWBA or adiabatic calculations. 

Finally the data in Fig. [11] for p + 17 O -> n + 17 F is 
analyzed in Ref. [H using DWBA. The DWBA calcula- 
tions leading to the ground state of 17 F are qualitatively 
similar to ours but fit the data better at the forward an- 
gles. For the transition to the l/2 + excited state of 17 F 
both calculations are quantitatively similar and fit the 
data equally well. 



V. CONCLUSIONS 

We have used the Faddeev/AGS three-body approach 
to study d + 12 C, d + 16 0, p + 13 C, and d + 17 reac- 
tions as a three-body system made up by a proton p, a 
neutron n and a structureless nuclear core A, the 12 C or 




0.1 , E p = 35MeV 
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FIG. 13: (Color online) Differential cross section for p + 13 C 
elastic and transfer reactions at E p = 35 MeV. Predictions of 
Model 2 with p-wave potential from Table [T] (solid curves) and 
IIIII (dashed curves) are compared. The experimental data are 
from Refs. jH [H, S3 . 



the 16 O. The interactions between pairs are the realis- 
tic interactions that describe np, nA, and pA scattering 
over the relevant energy range. The Coulomb interaction 
between the proton and the nuclear core is included in a 
numerically exact (converged) way. 

The aim of the present work is to demonstrate the 
possibilities and the shortcomings of the Faddeev/AGS 
three-body approach that provides, simultaneously, pre- 
dictions for all possible reactions, i.e., elastic, transfer 
and charge exchange such as, for example, p + 17 — > 
p + 17 0, p + 17 -> d + 16 0, and p + 17 -> n + 17 F, 
or d + 16 -> d + 16 0, d + 16 -> p + 17 0, and 
d+ 16 -> n + 17 F. 

Three different models (Ml, M2, M3) are studied in- 
volving energy-independent and energy-dependent opti- 
cal potentials that fit the nA and pA elastic scattering 
and whose parameters are fixed at a chosen energy or are 
allowed to vary over the energy range of the interacting 
pair, respectively. In the case of energy-dependent opti- 
cal potentials these become real at negative energies and 
support a number of single particle states that character- 
ize the (An) or the (Ap) nucleus. 

The results of our calculations indicate that transfer 
and charge exchange reactions at small angles are rather 
insensitive to the chosen model, but the elastic scatter- 
ing cross sections are highly sensitive to the choice of en- 
ergy dependence of the optical interaction (Ml versus M2 
and M3). Comparison with published CDCC, DWBA, 
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CCBA, and adiabatic calculations indicates that these 
approximate methods provide, in general, a better fit of 
the data than our calculations but are qualitatively sim- 
ilar to our results, particularly the ones of the "hybrid" 
model M3 that uses a partial- wave dependent optical po- 
tential whose parameters are energy-independent except 
in the partial waves that support the single particle states 
of the (An) and (Ap) nuclei. 
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APPENDIX A 

We present here selected results obtained with an alter- 
native N- 12 C p-wave potential whose spin-orbit strength 



is the same as in the other partial waves (see Table HJ . 
The strength of the central part is adjusted to reproduce 
13 C and 13 N lpi/2 ground state energies as in Table HT1 
thereby resulting different binding energies for the Pauli 
forbidden lp^n bound states that are given in Table HTT1 
together with the new values of potential parameters. 
The predictions of Model 2 with p-wave potential from 



v R (nA) v R (pA) 


V so (NA) 


lP3/ 2 ( 13 C) lp 3 / 2 ( 13 N) 


49.61 49.11 


5.5 


8.587* 5.507* 



TABLE III: Parameters of the alternative N-C p-wave po- 
tential together with the resulting binding energies for the 
Pauli forbidden bound state. See Tables HI and ILT1 for 

further explanations. 

Table HI and IIIII are compared in Fig. [13] for p + 13 C elas- 
tic and transfer reactions. The differences are rather in- 
significant when compared to the discrepancies between 
theory and data and therefore do not change the conclu- 
sions of this paper. Differences of similar magnitude can 
be seen also for the observables of d + 12 C reactions. 
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